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Abstract 



A qualitative study of the lightest glueball states in Euclidean SU (2) Yang-Mills theory quantized 
in the maximal Abelian gauge is presented. The analysis is done by generalizing to the maximal 
Abelian gauge the so-called replica model, already successfully introduced in the Landau gauge. As it 
will be shown, the gluon and ghost propagators obtained from the replica model are of the same type 
of those already introduced in fll, whose behavior turns out to be in agreement with that available 
from the lattice data on the maximal Abelian gauge. The model turns out to be renormalizable to all 
orders, while enabling us to introduce gauge invariant composite operators for the study of the lightest 
glueballs J^'-^ = O"'"^ , 2++ and ^ . The spectral representation for the correlation functions of these 
operators are evaluated to the first order, and the corresponding spectral densities are shown to be 
positive. Under the assumption of Abelian dominance, it turns out that the hierarchy for the masses 
of the lightest glueballs in the maximal Abelian gauge is in agreement with that already obtained in 
the Landau gauge, a feature which provides evidence for the gauge independence of the spectrum of 
the theory. 
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1 Introduction 



The maximal Abelian gauge |2l [3l H] is extensively employed in the study of nonperturbative aspects of 
Yang-Mills theories. This is the case, for example, of the dual superconductivity mechanism for confine- 
ment O ini [2] , according to which the low energy region of Yang-Mills theories should be described by an 
effective Abelian theory in the presence of monopoles. A dual Meissner effect arising as a consequence of 
the condensation of these magnetic charges might give rise to the formation of flux tubes which confine 
quarks. 

The possibility of achieving a description of the low energy regime of Yang-Mills theories in terms of 
an effective Abelian theory, also referred to as the Abelian dominance hypothesis [3 [U |9] , is being inves- 
tigated since several years, from both theoretical and lattice simulations points of view, see also \10\ lllj) 
for more recent investigations. In the maximal Abelian gauge, the Abelian degrees of freedom are iden- 
tified with the diagonal components of the gauge field corresponding to the diagonal generators of the 
gauge group. The remaining off-diagonal components decouple at low energies due to a dynamical gener- 
ation of an off-diagonal mass, a feature which has already received analytical |12^ \TE[ E] and numerical 
evidence [El m HTj . 

In this work we pursue the investigation of the nonperturbative aspects of Yang-Mills theories in the 
maximal Abelian gauge by addressing the issue of the glueball spectrum and of its comparison with 
the results already obtained in the case of the Landau gauge |18l ITU] . This goal will be achieved by 
combining several ingredients. Firstly, we shall generalize to the case of the maximal Abelian gauge the 
so called replica model, already introduced in the Landau gauge [20J. As we shall see, the replica model 
exhibits two nice features which are particularly useful in order to analyse the spectrum of the lightest 
glueballs. The first property concerns the gluon and ghost propagators, which turn out to be of the same 
kind of those already introduced in [1] and which are in agreement with the available lattice data [T7] . 
Moreover, unlike the model introduced in [1], the replica model enables us to make use of genuine local 
gauge invariant composite operators whose correlation functions can be employed to study the glueball 
spectrum. In addition, and in complete analogy with the case of the Landau gauge, the replica model of 
the maximal Abelian gauge will be proven to be renormalizable to all orders. As such, the model displays 
useful nonperturbative features. Secondly, we shall establish a suitable mass formula for the glueballs 
by relying on a kind of inspired SVZ sum rules, already employed in the case of the Landau gauge [18j . 
The combination of these two ingredients will enable us to provide a qualitative analysis of the location 
of the masses of the lightest glueball states. 

Having thus at our disposal a model which gives rise to diagonal and off diagonal gluon propagators 
in agreement with those observed in lattice simulations, we employ it in order to study the correla- 
tion functions of the local composite gauge invariant operators with the quantum numbers of the three 
lightest glueballs, namely J^*^ = 0'*''^,0 '",2++. In particular, thanks to the analytic structure of the 
diagonal gluon propagator, the correlation functions of the composite operators corresponding to the 
states J^'^ = 0++,0-+,2++, can be cast in the form of a Kallen-Lehmann spectral representation with 
positive spectral densities, enabling us to provide a qualitative study of the mass spectrum of the lightest 
glueball states in the maximal Abelian gauge. It turns out that, under the assumption of the Abelian 
dominance [3 El |9] , the resulting glueball spectrum is in agreement with that already obtained in the 
case of the Landau gauge [IBl [l9]. This result can be taken as evidence of the nonperturbative gauge 
independence of the physical spectrum of the theory. 

The paper is organized as follows. In Sect. 2 we give a brief account of the available lattice numerical 
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simulations on the gluon propagator in the maximal Abelian gauge. In Sect. 3 we construct the replica 
model for the maximal Abelian gauge and we show that the resulting gluon and ghost propagators have 
the same qualitative behavior of those reported in lattice simulations. In Sect. 4 we introduce the gauge in- 
variant composite operators corresponding to the three lightest glueball states J^*^ = O^"*" , , 2^~^ , and 
we determine the Kallen-Lehmann spectral representation for the corresponding correlation functions. 
Sect. 5 is devoted to present our results about the mass spectrum of the states J^*^ = 0+^,0 '",2++. 
In Sect. 6 we collect our conclusion. All technical details about the proof of the renormalizability of the 
replica model in the maximal Abelian gauge can be found in the Appendix A, while in the Appendix B 
the plots of the mass hierarchy of the states J^'-^ = 0++, ^, 2++ are displayed. 



2 A look at the lattice data on the gluon propagator in the maximal 
Abelian gauge 

Before giving a short account of the lattice investigations, let us remind that, in the SU{2) maximal 
Abelian gauge, the gauge field A/j, is decomposed into diagonal and off-diagonal components, according 
to 




where ((j",(T'^) stand for the Pauli matrices. The diagonal component, corresponding to the diagonal 
generator cr^, is A^, while a = 1,2, denote the off-diagonal components. In the maximal Abelian 
gauge, the fields {A^,A^) are constrained by the following conditions |2l|3lll] 

DfA'^ = d,A;-ge'''A,A'^ = 0, 
d^A^ = 0. (2) 

It is useful to remind here that the non-linear maximal Abelian gauge condition, D^A^^ = 0, can be 
implemented by minimizing the auxiliary functional [Ij 

F[A] = I d'xA^Al , (3) 

with respect to the gauge transformations. The functional ([3]) possesses a lattice version [161 117]., mak- 
ing possible the numerical implementation of the maximal Abelian gauge by selecting the configurations 
which correspond to local minima of the auxiliary functional. 

The infrared behavior of the gluon propagator in the SU{2) maximal Abelian gauge has been investi- 
gated by several authors through lattice numerical simulations. In [15], both the diagonal and off-diagonal 
propagators were analysed in configuration space, exhibiting a behavior in agreement with the Abelian 
dominance. The off-diagonal gluon propagator turned out in fact to be suppressed with respect to the 
diagonal propagator. Later on, in [16j, the diagonal and off-diagonal gluon propagators have been studied 
in momentum space. Also here, the off-diagonal propagator is suppressed with respect to the diagonal one. 

More recently, this issue has been addressed by the authors |17j on a larger lattice. Both gluon and 
ghost propagators have been analysed, and particular attention has been devoted to their behavior at 
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very low momenta p ^ 0. According to ^17] . the transverse diagonal propagator can be fitted by 



Dd^aaip') = 4 > (4) 

a + bp^ + cp^ 

were the parameter rridiag = s/j, is found to be rridiag ~ 0.72GeV. It is worth to point out that the 
diagonal form factor D^iagiP^) does not vanish at zero momentum, p = 0, and that expression Q is 
precisely of the same kind of that found in the case of the Landau gauge within the so called Refined 
Gribov-Zwanziger framework [2H [22l [23] as well as within the Landau replica model [20]. It is worth 
underlining that the nonvanishing of the gluon propagator at zero momentum in the Landau gauge has 
been also seen in lattice simulations, see for example [2^[25 l [M l l27j. As we shall see in the following, this 
is a very welcome feature which will play an important role in the discussion of the glueball spectrum in 
the maximal Abelian gauge. 

The best fit for the transverse off-diagonal gluon propagator is of the Yukawa type |T7] , namely 

where the parameter rrioff = takes now the value ?ti.o// =~ 0.97Gey. Although, as expected from 
the hypothesis of Abelian dominance, the value of the off-diagonal parameter rriofj is greater than nidiag, 
a clear evidence for the Abelian dominance would demand a more complete analysis, perhaps using larger 
lattices. 

Although not needed for the purposes of the present investigation, it is worth to mention that the 
ghost propagator decomposes into symmetric and anti-symmetric components in color space, i.e. 

{(f{p)c\-p)) = S'^'^Dsymm + e^^Dasymm • (6) 



It is interesting to observe that the best fit for the symmetric form factor Dsymm obtained in [17] is also 
of the Refined Gribov-Zwanziger type for the maximal Abelian gauge, as obtained in [Ij, namely 



1 + dp^ 
a + hp^ + cp^ ' 



with a = 0.45(l)Gey2, b = Ll(3), c = 0.73(30)Gey-2, d = 2.l{9)GeV~'^ . In particular, there is no sign of 
infrared enhancement in expression ([T]), which attains a finite non- vanishing value at p 0. Concerning 
the antisymmetric component, Dasymm, the authors [T7] have provided a first preliminary evidence that 
it might be non-vanishing, a fact that could be linked to the existence of certain non-vanishing ghost 
condensates of dimension two which can be naturally introduced in the maximal Abelian gauge, as argued 
in [28]. For specific details on the potential errors related to finite size and discretization effects we remind 
to the original lattice references [El [Ml [IT] . 
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3 Construction of the replica model in the maximal Abelian gauge 



In order to construct the replica model, let us start with the Faddeev-Popov action for SU (2) Yang-Mills 
quantized in the maximal Abelian gauge, namely 



+ibd^,A^ + cd^id^c + ^e^'^uc')! , (8) 



S[A] 

where (-F^^,, F^y) are the off-diagonal and diagonal components of the field strength 

f;, = DfAl - DfA'^ , F,, = d,A, - d,A, + ge'^^AlAl , (9) 
and D'^f' is the covariant derivative respect to the diagonal component, eq.Q, 

Df = 6^'d^-ge-'A^. (10) 

The indices a, b, c appearing in eqs.Q, ([9]), (10) refer to the off-diagonal components, and run from 1 to 
2, i.e. a,b,c = 1,2. The fields {b^,b) are the Lagrange multipliers implementing the gauge conditions 
([2]), and (c",c"), (c, c) are the off-diagonal and diagonal Faddeev-Popov ghosts. Also, denotes the 
Faddeev-Popov operator of the maximal Abelian gauge, being given by 

^ab ^ -D^^^Df - g^e^^e^'^A'l^A'l . (11) 

The action ([s]) is left invariant by the nilpotent BRST transformations: 

sA^ = -{Dfc'' + ge'^'^A^^c) , sA^, = -(9^c + ge^^Alc^) , 
sc"" = ge'^^c^c , sc = |e"*c"c'' , 

sc" = , sc = ib , 
sb'' = 0, sb = 0, (12) 

with 

sS[A] = . (13) 

Furthermore, it is known that, for renormalization purposes \29\ I30j . an additional BRST exact term 
proportional to a gauge parameter a has to be added to the action S[A], which is replaced by 

So[A] = S[A] + , (14) 

where 

'a 
2 



J (fx {b^b" + 2ige''^ b^'c^c + g^(f(^>c''c^^ . (15) 



The action S'of^] enjoys multiplicative renormalizability to all orders [3D]. The limit a — )• has to be 
taken after the removal of the ultraviolet divergences. 
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We are now ready to start the construction of the rephca model in the maximal Abelian gauge. To 
that end, we shall follow the procedure already outlined in the case of the Landau gauge [201 [IB]- We 
first introduce a replica of the action S'o[^], namely 



So[B] 



(fx 



1 "(a r<a \ ^ r< r< i aUO. nabTyb -a k a'^^ , .b 



(16) 



with 



M 



■ab 



DfBl-DfB'^, 

d^B, - d,B,, + ge'^^'B^^Bt , 

-Dl'Df-g\'^'e'''BlBl. 



(17) 



The fields (B^^, B'^,uj,uj,uj"',uj°',b,b°') are a replica of those of expression (14). The action So[B] is left 
invariant by the BRST symmetry 



sB'^ = -{Dfu^ + ge'^^Bf^io 



-ab d6 , 



sB„ = -{d^uj + ge^'>B^^uj'>), 



SLO 



ge^^uj^uj , 



= , 



SUJ 



2 



SUJ = ib , 
sb = 0, 



(18) 

sSo[B] = . (19) 

Further, following the case of the Landau gauge [201 [18], the two expressions So[A] and So[B] are coupled 
through a soft term 

S^ = iV2^^ j d^xA^B^ , (20) 

where i? stands for a a free mass parameter. Therefore, we consider the action 

S = So[A] + So[B] + . (21) 

It is worth to point out that 5* enjoys the following discrete symmetry which, for obvious reasons, will 
be called the mirror symmetry: 

{A';,,A,) ^ {b;,b,), 

(c^c) o (a;^w), 

(6^6) o (6^6). (22) 

Moreover, due to the presence of the term S^, expression ( |21| ) is not left invariant by the BRST trans- 
formations which, as in the case of the Landau gauge [20[ I18j . turn out to be softly broken, i.e. 



sS = sS^ = -iV2^^ / d^x \{d^c + ge''''A''^c'')B^, + {d^,uj + ge^''B^^oj^)A^ 



(23) 
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Nevertheless, as we shall see in the following, the presence of the BRST soft breaking does not jeopar- 
dize the renormalizability of the model as well as the introduction of local BRST invariant composite 
operators describing the lightest glueball states. As already argued in [20], the mechanism of the BRST 
soft breaking could be suitable in order to achieve a good description of the infrared region of a confining 
gauge theory. On one hand, the soft parameter -i?'^ induces an infrared modification of the propagators of 
the elementary gauge fields in such a way that they cannot describe the propagation of physical particles. 
On the other hand, a set of suitable composite BRST invariant operators can be introduced in such a 
way that their correlation functions can be employed to study the physical spectrum of the theory. This 
framework is precisely what one could expect in a confining theory like pure Yang-Mills theory where 
gluons are not part of the physical spectrum, which is given by glueballs. 



Before addressing the issue of the introduction of the glueball operators, it is worth to elaborate more 
on the properties of the BRST soft breaking. In particular, in the next section we shall prove that the 



breaking ( 20 ) can be cast in the form of a linear breaking, i. e. of a breaking which is purely linear in the 



quantum fields. This is a relevant property with far reaching consequences. It is known in fact that a 
linearly broken symmetry can be directly translated into a powerful Ward identity |31] , which will enable 
us to give a purely algebraic proof of the renormalizability of the model, as recently done in the case of 
the Gribov-Zwanziger theory, see \32\ I33|. 



3.1 Converting the soft breaking into a linear one 



In order to convert the soft breaking, eq.(23), into a linear one, we follow the procedure outlined in 



[32| [33] and introduce a BRST quartet of auxiliary fields 



sX = a . 
sa = X. 



sa = , 
sX = 0. 



(24) 



Next, we consider the following action: 

Sa = s j (fx [xA^B^ - Act + ^Act^ + 1^/2^"^ j d^xa{a 



j (fx {aA^B^ + ^ct2 - ct(ct - iV2i?2) 



+A 



A + {d^c + g£'''Ay)B^ + {d^uj + ge^''By)A 



ab Tja, ,b\ 



(25) 



where the term proportional to the dimensionless constant parameter is allowed by power-counting. 
It is easy to check now that the action So- breaks the BRST symmetry in a linear way, i.e. the resulting 
breaking term turns out to be linear in the fields, namely 



sSo = i\/2??2 j (fx A(x) . 



(26) 



It remains now to prove that S^ is equivalent to S^. To that end, let us notice that the auxiliary field 
ct(x) plays the role of a Lagrange multiplier enforcing the constraint 



CT - iv/2??2 = . 



(27) 
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In fact, integrating out the field o", we have: 



(28) 



We can now perform a change of variables (A, A): 

A' = A, 



A' = ~X + id^c + ge^'A';,c'')B^ + {d^u; + ge''''By)A 



fj. ■ 



with unity Jacobian. Therefore 



S^ = J (fx (^^^^ A^B^ -Ci^^ + A'A' j . 



(29) 



(30) 



As the last term, A' A', is completely decoupled from the rest of theory, we can neglect it. Up to a constant 
vacuum term, Ci '&^, we have shown that is equivalent to S^. Thus, the quadratic soft breaking (23) 



has been converted into a linear breaking. From now on, we shall consider the following equivalent action 



Sun = So[A] + S^[B] + S^. 



(31) 



Let us end this section by remarking that the equivalence between the two formulations, i.e. the softly 
and the linearly broken formulation, can be promoted at the quantum level. In fact, by following the 
same steps outlined in the case of the Landau gauge [32l |33] , it is not difficult to establish the formal 
equivalence between the partition functions of the two formulations. 



3.2 Introducing dimension two operators 



We can pursue now the task of reproducing the infrared behavior of the gluon and ghost propagators 
observed in lattice simulations, eqs.Q, ([5]), ([T]). The action (31) is not yet in its final form in order to 
fulfill this task. In fact, we can introduce a whole set of dimension two operators, 0°^J , 0'^2^ , Occ, Ocxc, 
which will enable us to reproduce in a nice way the lattice results, while giving rise to an all order 
renormalizable action. These operators read: 



off 
A2 



diag 
Or., 



Or 



^{A^AI + BIBD, 
]^{A^A^ + B^B^), 

ge''\(fc^ + oj^J) . 



(32) 



As done in the case of the soft term of eq.(20), all operators (32) can be introduced in such a way 



off 



that the resulting breaking of the BRST symmetry is linear. Consider, for example, the operator 
Proceeding as in the previous section, we introduce a BRST quartet 



sip = p, 
sp = (p, 



sp 

s(p 



0, 
0, 



(33) 
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and write the following term 

Sp = s J (fx 0°lJ — 'PP+ y V'P^ + y d^x p{x) 

= J d'x (pO%^ + |/ - p(p -m') + ^(^^- sO%f) ) , (34) 

where m is a mass parameter and ^2 a dimensionless constant allowed by power counting. Expression Sp 
gives rise to a linear breaking of the BRST symmetry, namely 

sSp = vr? J d^x (p{x) . (35) 

Moreover, Sp is equivalent to the introduction of the term / d^x 0°JJ . In fact, up to a constant term, 
we have 



.2 /no// , C2^ 



Sp = j d^x\rn'0"l^ ^^-^m"^ , (36) 

where use has been made of the constraint p = rr? and of the change of variable (p ^ (p — sC^J" , which 
has unity Jacobian. The remaining operators can be introduced in the same way. Denoting by {tp, h, h, tj)) 
the BRST quartet needed for the introduction of the operator Occ, we have 

st/j = h , sh = , 

sh = '4), sijj = 0, (37) 
so that the corresponding term accounting for the introduction of Occ is given by 

Sh = s j d^x Occ -i}h+ ^'^hj J d'^x h{x) 

= j d\ (^hOcc + - h{h -M^)+^{'tjj- sOcc) ) , 
= Jd'xI^M'O-cc+^fM^y 
Analogously, for the operator O^j"^) we make use of the quartet 



(38) 



sx = V, ST] = 0, 

sf] = X, sx = 0, (39) 

and write the following term 



Sn = s j d^x ^X^^l^^ -X'7+ yX^J ^ J d^^^vix) 

= J d'x [rj Of^^ + |r?2 - fj{ri -p')+x{x- sOf^^) 
= j d^xlp^O"'^' + '^-^A. (40) 



Finally, for Ocxc; we have 



su = ^ , = , 

sC = u, su = 0, (41) 
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and 



uh + ^ui ] +v'^ I d^x^ix) 



d^x ( e Oexc + - e(e -v^) + u{u- so, 

d''x{v'0,x, + ^vA. 



(42) 



3.3 The final expression of the classical action for the replica model in the maximal 
Abelian gauge 

We have now all ingredients to write down the complete classical action describing the replica model in 
the maximal Abelian gauge. First, we collect all terms introduced in the previous sections and consider 

Sun + Sp + Sh + + . (43) 

Moreover, in order to discuss the renormalizability of the model, we need to introduce external sources 
coupled to the composite operators corresponding to the nonlinear BRST transformations [3lj, namely 

Ssources = J d^x [Q"^ (sA"^) + Cl"^ {sB"^) + (sAf,) + Cl^, {sB^) + {sc") + (SUJ^) + L (sc) + L (suj)] . 

(44) 

where (O^, Q'^, f^^, O^, L°, L", L, L) are BRST invariant external fields. It is apparent that the mirror 



symmetry, eq.(22), can be immediately extended to the external sources: 



{L^L) o (L^L). (45) 
Finally, it turns out that the following term 

Sextra = j d'^x {ki (7 p + k2 cn] + ah + k4pri + k^ ph + /cg rjh) , (46) 

is required by power-counting. The parameters {ki, k2, k^, k^, k^, k^) will enable us to correctly take into 
account the possible mixing arising at the quantum level among the various operators of dimensions two 
which have been introduced, as it will be discusses in Appendix A. 

Therefore, the final form of the classical action describing the replica model in the maximal Abelian 
gauge is _ 

5^ = Slin + Sp + Sfi + S'q + S"^ + Ssources + Sextra ■ ('^''') 

It exhibits the important feature of breaking the BRST symmetry in a linear way, i.e. 



d^'x (iV2^'^ X{x) + (p{x) + i){x) + p^ x{x) + u{x)^ . (48) 



As we shall see in details in Appendix A, eq.(48) can be directly translated into a system of helpful 
Slavnov- Taylor identities. Moreover, in addition of the Slavnov- Taylor identities, the action E exhibits a 
huge set off additional Ward identities which will ensure the renormalizability to all orders. 
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3.4 The gluon and ghost propagators from the rephca model and their comparison 
with the lattice data 



Having identified the final form of the action of the rephca model in the maximal Abelian gauge, eq.(47), 
let us proceed by evaluating the gluon and ghost propagators and see how they compare with the corre- 
sponding lattice expressions, eqs. Q, ([s]), ([T]). 

For the transverse off-diagonal gluon component we obtain a Yukawa type propagator 

{A';,{k)At{-k)) = {Bl{k)Bl{-k)) = - ^) ^"^^ • (49) 

The diagonal gluon propagator turns out to be of the type obtained in [1]. This kind of propagator also 
appears in the Refined Gribov-Zwanziger theory of the Landau gauge, see [21 ^ 122 1 [23]. i-^- 

{A,{k)M-k)) = {B,{k)B^{-k)) = - ^) • (50) 

For the mixed A — B propagators we get 

and 

{A-ik)Bt{-k))=0. (52) 
The symmetric off-diagonal ghost propagator is given by 

{C^{k)c\-k))syn.m = mk)u'{-k))symm = + ^2)2 + ^4 ' (^3) 

while for the anti-symmetric off-diagonal ghost we obtain 

a b a b ab 

(c {k)c { — k))anti—symm — {'-^ {k)LO { — k)) anti— symm — j!j2 _j_ _/y2^2 _j_ ^4 ^ ' (^4) 



Looking at expressions (49), (50), (53) we see that they are precisely of the same qualitative kind of 
those reported in lattice simulations, and already obtained in [1]. This feature shows that the replica 
model might be very helpful in order to investigate the infrared properties of confining Yang-Mills theories. 



Before ending this section, it is worth to spend a few words on the massive parameters (t?"^, m"^, M'^,ij,^, v"^^ 



which appear in the classical action (47). In the present work, these parameters are regarded as free 
parameters whose introduction can be supported by the requirement of agreement with the behavior of 
the gluon and ghost propagators obtained in lattice simulations. Moreover, although being out of the aim 
of this work, we point out that a possible framework to provide a deeper meaning to those parameters 
could be achieved by implementing, within the framework of the replica model, the restriction of the 
domain of integration in the path integral to the so-called Gribov region Q, of the maximal Abelian 
gaug^ In fact, as it stands, the action of the replica model, eq.(47), is still plagued by the presence 



^We remind here that the Gribov region Q, in the maximal Abelian gauge can be defined as the set of all configurations 
for which the Faddeev-Popov operator of the maximal Abelian gauge is strictly positive, namely 

^ = {A^,Al, a^A^ = 0, DfAl^Q, 0}. (55) 



11 



of Gribov copies. This observation follows by realizing that the replica model has been constructed 
by making use of the maximal Abelian gauge conditions Wj^A^^ = and D^B'^ = 0. Consequently, 

the Faddeev- Popov operators ^A"'^ and ^A"''' in eqs.(ll),(17) are nothing but the usual Faddeev-Popov 
operators of the maximal Abelian gauge. As such, they have zero modes which correspond to Gribov 
copies. A way to take into account the existence of these copies would be that of implementing, in the 
case of the replica model, the restriction in the domain of integration to the Gribov region of the maximal 
Abelian gauge, as discussed in EHl EH EH]) see also [2D] for a treatment of the Gribov issue in the 
replica model for the Landau gauge. This would give rise to a set of gap equations which would provide 
a dynamical origin for the parameters (i?^, m^, M^, u ^). Though, for the time being, these parameters 
will be treated as free parameters. 



4 Composite operators and spectral functions 

As already pointed out in the introduction, one relevant feature of the replica model is the possibility 
of introducing BRST invariant composite operators whose two-point correlation functions exhibit the 
Kallen-Lehmann spectral representation with positive spectral densities. As such, these correlation func- 
tions can be given a physical meaning, being suitable for the investigation of the spectrum of the glueballs. 

We underline that the introduction of these BRST invariant composite operators is a nontrivial fea- 
ture of the replica model, given that these correlation functions have to be evaluated with confining 
propagators which have complex poles as, for example 

It is worth mentioning here that this type of propagator also exhibits positivity violation, see [22) . 



Let us also remind that the introduction of BRST invariant composite operators exhibiting good spec- 
tral properties is still an open problem of both Gribov-Zwanziger |391 HQ] and Refined Gribov-Zwanziger 
[m [221 [23] frameworks. 

The mechanism which enables us to construct composite operators which, despite the use of a con- 
fining propagator which has complex poles, give rise to correlation functions with nice spectral properties 



has been elucidated in [JT] and relies on the concept of i-particles. From expression (56), it is apparent 
that the confining gluon propagator in the right hand side can be seen as describing the propagation of 
two unphysical modes with complex conjugate masses (^^ it i\f2'd'^). These unphysical modes have been 
called i-particles in |i41j . Moreover, it turns out that, for a certain type of composite operators, only pairs 
of i-particles with conjugate masses contribute to the two-point correlation functions of these operators. 
In this case, the complex conjugate poles combine each other in such a way that a cut along the real 
axis in the complex fc^-plane emerges at the end of the computation ^T|, allowing us to write down a 
nice spectral representation. In the case of the replica model, it turns out that this particular class of 
composite operators is precisely given by the BRST invariant operators. This is what one would expect 
in a confining theory: the elementary fields do not correspond to excitations of the physical spectrum, 
described by a suitable set of BRST invariant composite operators corresponding to bound states of the 
confined unphysical elementary excitations. 

Let us proceed thus by introducing the BRST invariant composite operators with the quantum num- 
bers of the three lightest glueballs. These states are classified according to the values of J^*^, where J is 
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the angular momentum, P the parity, and C the charge conjugation, see |42j for a recent review on glue- 
baUs. The hghtest glueball states have J^*^ = 0^'*', ^, 2^+. Following the same procedure of the replica 
model of the Landau gauge [IB], for the composite operators corresponding to J^'-' = ^, 2~^~^, we 
have 



'pA pA _qAqA\ 



1 

2 

— ( TP"- TP°- — r""- C"- J- TP TP — C C \ 



(57) 
(58) 



Oo-+(x) 



[pApA _ qA qA ^ 



i^uy^pa + Pp.vFpa 



t^pvpa\^ p^yJ- pa 



^ \JLV^ pa 



[C'2++(x)]^j^ - {^paPv^ - l^PpLvPa^ [F^^Pf^ - G^^G^^ 



PpaPufS - ^PpuPaP ) {FaaPpa + PaaFpa - G^^G^^ - GaaGfSa) 



(59) 
(60) 

(61) 
(62) 



where Pp^ = Sp^d'^ — dpdy is the transverse projector. As pointed out in [331 US]) the last operator 



generates a pure 2"^''", as follows by noticing that expression (62) is symmetric, traceless and conserved. 



It remains now to evaluate the two-point correlation functions 

{Oi{k)Oi{-k)) , i = 0++,0-+,2++ , (63) 
and show that they can be cast in the form of a Kallen-Lehmann spectral representation, i.e. 



(0,(k)Oi{-k)) 



Making use of ^-particles, see eq.(56), for the correlation functions (63) at one- loop order, we get 



{Oi{k)Oi{-k)) 



fi{p,k-p) 

(27r)4 (p2 + ^2 + i^^2^(p2 + ^2 _ ^^^2) 



(64) 



(65) 



where fi{p,k — p),i = 0^''',2++,0 + are polynomials in the scalar products of the momenta {k,p). There 



are several ways to cast expression (65) in the form of a Kallen-Lehmann representation. As recently 
discussed in [53], a powerful framework is that of employing Cutkosky's rules, and performing an analytic 



continuation to complex masses in Euclidean space. Noticing that, at one-loop order, the correlation 
functions {Oi{k)Oi{—k)) can be decomposed into diagonal and off-diagonal components, namely 



{Oi{k)0^{-k)) = {0f''\k)0f''\-k)) + {0f^ik)0ffi-k)) 



(66) 
(67) 



it turns out that, after a rather long computation, the corresponding spectral representations are given 
by the following expressions 



{O°ff{k)0°ff{-k))o++ 



dT 



1 



8vr2 74^2 r + A;^ 



1 



4m^ 



(4m^ + 2r^ + 8m^{m^ - t)) , 



(68) 
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1 



2++ 



dT 



Gtt^ 74^2 r + /c^ V T 



1 - — 2T\6m^ - SmV + r^) , 



(69) 



1 



dT 



1 (r^ -4m^ + W(m2 -r)) , (70) 



{0''^'{k)0'^'^^{-k))o.. = ^ r , -4^a/i - ^ - ^(V + + 2r' + - r)) , 

(71) 



Stt^ 72(/.2+VM4+2i,4) T + A;^ V r 

+3r2(6//^-3^V + T2) + 32i?*) , (72) 



(0'^^"S(A;)0'^*'^f(-A:))o- 



(73) 



It is worth observing that all spectral densities entering expressions (68)-(73) are positive definite within 
the corresponding range of integration. 



5 Qualitative analysis of the spectrum of the lightest glueballs, J = 

0++,2++,0-+ 

With the information encoded in the Kallen-Lehmann spectral representation of the correlation func- 
tions of the glueball operators, we can provide a qualitative study of the spectrum of the lightest states 
J^^ = O"*""*", 2"'""'", in the maximal Abelian gauge. As we shall see later on, the hypothesis of Abelian 
dominance will play a pivotal role. 

The first task to study the spectrum of the glueballs is that of establishing a suitable mass formula 
through the use of the spectral densities evaluated in the previous sections. To that end we shall adopt 
the set up which has been successfully developed in the case of the Landau gauge [T8|, based on a kind 
of SVZ-type sum rules [ISl 061 HTj. For the benefit of the reader, let us give here a short account of the 
method. As it is customary in the SVZ approach to QCD [38], we start by considering the two point 
correlation functions 

n,(g2) = j d^xe'^^{0,{x)Om) , (74) 

where 0%, i = 0'^~'',2^^,0 stand for the local composite gauge invariant operators which generate 
glueball states with quantum numbers J^^ = 0'*'"'' , 2"*"^ , ^. 
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On physical grounds, a truly nonperturbative evaluation of Hi^q^) would enable us to write an exact 
Kallen-Lehmann spectral representation 

Ui{q^) = - dr 75 

which is expected to follow from the unitarity and analyticity properties of the underlying nonperturbative 
theorjj^ We thus proceed by employing a one-resonance parametrization for /mnj(r), namely 

^ = n, 6{t - m\) + Q(t - tI)p,{t) , (76) 

vr 

where mf denotes the glueball mass in the i-th channel, Tq is the threshold for the continuum part of the 
spectrum, and /?i(r) the corresponding positive spectral density. Of course, the real values of TZi, mj, 
Tq and of the spectral density /Oi(r) are unknown. So far, the best which can be done is computing the 



correlation functions (74) by trying to encode as much nonperturbative effects as possible, as summarized 
by the following equation 



..„.2+ , "-7X^=1-, (77) 



where 11 ■ ^ stands for the expression of the correlation function ( 74 ) which one has been able to evaluate 



in practice. Expression (77) establishes the so-called sum rules, enabling us to give estimates of the 



glueball masses in terms of the nonperturbative parameters present in H"^. As done in [18], we shall 



attempt at evaluating the right hand side of eq.(77) by employing a suitable trial action Striai which 



encodes nonperturbative information about the infrared dynamics of Yang-Mills theories. Of course, the 
choice of Striai is a nontrivial matter which has to account for the following requirements: 

• i) the action Stnai has to display a nonperturbative character, encoded in the presence of a set 
of nonperturbative parameters {A}. Moreover, in the limit A = 0, in which these parameters are 
removed, the action StriaiW has to reduce to the usual perturbative Faddeev-Popov action Spp, 
i.e. 

Striai W _ = Spp , (78) 

• a) it accounts for gluon confinement. This means that the two point correlation function of the 
elementary gluon field evaluated with Striai W cannot be cast in the form of a spectral representation 
with positive spectral density, so that it cannot be interpreted as the propagator of a physical 
particle. 

• iii) Striaii^) has to be renormalizable, meaning that consistent perturbative calculations can be 
worked out. 

• iv) StriaiW should enable us to introduce gauge invariant or, equivalently, BRST invariant local 
composite operators {Oi} with the quantum numbers J^'^ = 0+"'" , 2"'""'" , '", whose two point 
correlation functions 

n7(g2) ^ nf-'(g2, A) = I d^xe^«-'(O.(x)O,(0)) , (79) 



We remind here that, in some cases, the spectral representation, eq.(75l, might require appropriate subtraction terms 
in order to ensure convergence. These terms are not written down, as they will be removed once the Borel transformation 
wiU be taken HSl. 
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exhibit the Kahen-Lehmann spectral representation, at least at their lowest order, i.e. 

retrial 



^7(g^) = ^r-'((z^A) 



(i)trial / 1 \ 

i-^i 



dr 



Pi 



T + 



(80) 



Requirements i)-iv) look quite stringent. As such, they might provide a satisfactory set up in order to 
achieve a useful predictive expression for the correlation functions 11 ■"'^'(g^. A), eq.(79). 



As trial action we shall adopt the action of the replica model of the maximal Abelian gauge given 
in the expression (47), i.e. 

StrialW = ^ , (81) 
so that the quantity n*'~*°'(g^. A) in eq.(80) will be identified with the corresponding expression evaluated 
with the replica action S, namely 



n: 



trial 



n: 



repl 



trial 



Pi 



P. 



repl 



(82) 

As shown in the previous sections, the action S fulfills all requirements i) — iv). In particular, the 
nonperturbative character of S is encoded in the mass parameters (i?^, m^, Af^, u^) which enable us 
to reproduce the infrared behavior of the gluon and ghost propagators observed in lattice simulations. 
Furthermore, taking into account that, according to eqs. (68)-(73), the spectral functions decompose 
at one-loop order as 

^repl _ ,^repl, , [p^^^^ f , (83) 



Pi 



it follows that eq.(80) becomes 



VPi \diag 



_(i)repZn ^ _|_ 

\diag ^ 



[pr'ir)]off 



{i)repl-i 

b loff 



dT 



(84) 



Before going any further and, as done in the SVZ sum rules 
estimate for the quantity 

Pi{t) 



dT 



HZl HH] , we need to provide an 

(85) 



appearing in eqs. (76), (77). Requirement iv) above states that the expression for II"^ in eq.(77) reads, at 
the lowest order: 



q2 _|_ 



+ 



T + q^ 



(i)repZ 1 
\diag 



30 r repl / \ i 

dr + 0{h^) . (86) 



Following [18j, we shall estimate expression (85) by assuming that, up to the required order, the exact 
spectral functions Pi(r) are well approximated by those which have been evaluated explicitly at one- loop 
order by means of the replica action S, that is 



P^{r) = [pr\r)Ua, e{r - a [r^'^-^'],,,,) + [p^ (r)]o// 9{r - h [r^'^--%;^) , (87) 
where a,b & stand for free parameters accounting for the difference between the value of the true 
unknown physical threshold Tq in eq.(85) and the thresholds ^[TQ*^^^^^]djagi [■'"o*^'''^^']o//^ which we have 
been able to evaluate in practice, eqs.(68)-(73). Therefore, for the final form of the sum rules which we 
shall employ, we write 



repl I 



{i)repl-i 



Hi 



q^ + mf 



{i)repl 1 
\diag 



+ 



{i)repl 1 

\ofl 



dT 



T + . 



This equation will be taken as the starting point in order to derive a mass formula for the glueballs and to 
provide a qualitative analysis of the lightest states J^^ = 0"'"+,2++,0 in the maximal Abelian gauge. 
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5.1 A mass formula for the glueballs 



All ingredients are now at our disposal to work out a mass formula for the glueball states. As in the SVZ 
sum rules approach [45t ITT] , we shall make use of the Borel transformation. Let us briefly illustrate 
how this works. One starts by considering an equation of the following type 



dT- 



0- T 



for some suitable F, rni and o"(r). Equation (89) is precisely of the kind of expression 



(89) 
Let us 



remind some basic properties of the Borel transformation B^, i.e. 



B 



+ TTl; 



(90) 



where /ig is the so called Borel mass, see for example liSl . Applying thus the operation to both sides 
of eq. ( 89 ) , we get 





_m2 














dr a{T 


)e 
























= 




dr T (J 


r) 













(91) 



so that one derives the mass formula [18 



P^dTT(j{T)e "e 

_ T 

J^^dTa{T)e ^ 



(92) 



Let us make use of this formula by focusing, without loss of generality, to the case of the glueball state 
0"*"^. Before applying the Borel transformation, it is useful to perform helpful changes of variables in the 



integrals of expressions (68) and (71). For the diagonal part, eq.(71), we shall set 



f = T-2fX% 



(93) 



so that 

(0'^^'^9(A;)0^*"9(-fc))o4 





poo 


87r2 


J A 




POO 


87r2 


J A 



df 
df 



1 



V7 



A2 



1 



;2f2 + 



fc2 + f + VA2 -8??V f + VA2 - 87?4 



Vf2 



A2 



f2f2 + A^ 



(94) 



For the off-diagonal component, eq.(68), we shall make use of 

f = T — 2m^ , 
r = 2m^ . 



(95) 
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Thus 



s^ir U' + ^ + ry f + r 



(2f2 + r2) . 



Applying now the Borel transformation to the diagonal component, one gets 



_{i+2i«f) 



A;2 



where we have introduced the following dimensionless quantities 



T + 2;u2) e = -A (t + ^) e-P^e-P" 



(96) 



(97) 



Analogously, for the off-diagonal part: 



where 



1 



fc2 

A;2 + f + r 



A 

P = ^ ^ 
2/^2 



(f+T) 
" ,,2 



(f + r)e "8 = -^(s + l)e-"e- 



(98) 



(99) 



^1 



Finally, putting all together, we obtain 

dt^/t^^ (2*2 + 1) e-P^e-P" + 34 ti dsV^^^ {2s^ + l) e'^'^e" 



0+^ 



{a,b,p,q,r) = A- 



(2t2 + 1) e-P*e-P^ + 3^ ds^^ {2s^ + 1) e— e" 



(100) 



(101) 



This formula follows by performing the following steps: we start from expression (88) which is a particular 



case of expression (89). In particular, the spectral densities appearing in eq.(88) can be read off from 



eqs.(94), (96). Further, the Borel transform in eq.(90) is applied to the specific case of the spectral 



spectral densities given in eqs.(94), (96), resulting in the mass formula (101). Similar expressions can be 
derived for the masses corresponding to the states 2"''"^, ^ . 



5.2 Taking into account the hypothesis of Abelian dominance 



Having been able to work out a mass formula for the glueball states, the next step is that of making a 
comparison among the three masses m2_^_,. , 171^++ , 'm'^_+ , and check out if their location is in agreement 
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with the results obtained in lattice simulations [49j . according to which the 0"*"^ is the lowest state, fol- 



lowed by the 2+"*", the being the heaviest one, i.e. 171^++ < m'^++ < rn^ 

At this stage, and before discussing the mass hierarchy, we need to take into account the consequences 
stemming from the hypothesis of Abelian dominance 13 [U [9], which is believed to be a key feature of the 
maximal Abelian gauge. Let us remind here that, according to the Abelian dominance, the off-diagonal 
degrees of freedom are expected to decouple at low energies, due to the dynamical generation of a large 
off-diagonal mass, a fact which has already received a certain amount of evidence, from numerical lattice 
simulations [T6l [T7] as well as from theoretical calculations \\.2\ [T3l . Even if the issue of the Abelian 
dominance is still an open subject, let us proceed by assuming that it holds and that, accordingly, the 
off-diagonal components acquire a large dynamical mass which decouple them at low energies. As we 
shall see, the hypothesis of the Abelian dominance will have great impact on the glueball mass formula 



( 101 ). Looking in fact at expression ( 101 ), one immediately recognizes that the contribution coming from 



the off-diagonal sector is encoded in the two integrals 



where 



^ / ds \/s2-l e-^' e-' , W ds — (2^2 + l) e"" e"'^ , (102) 

P Ji p-^ Ji s + 1 

r m? 

r = ^=2^, (103) 

and m stands for the mass of the off-diagonal components of the gauge field. If we assume now that, 
according to the Abelian dominance hypothesis, m is larger enough than the remaining mass parameters 



(/i, ??) which enter the expression of the diagonal gluon propagator, eq.(50), it follows that the off-diagonal 
contribution to exp 
g-r(s+i) appearing 

state 0"*"+ becomes 



contribution to expression (101) turns out to be very much suppressed, due to the exponential factor 
g-r(s+i) appearing in eq.(102). Neglecting thus the off-diagonal contribution, the mass formula for the 



mQ++^ Abelian dom. = A '—— -== . (104) 



Similarly, for the other states 2"'""'" and O"*"^, we obtain 

dt {q + t)2(5g2(2t2 + i) + 155^(^2 + 1) + gt^ + 8t + I) e~P' 

t+q 



mL+(a,p) = A ^^^^^^ ' ' ' ^"^^^^ ' ' - i ^ V- — ^ ^-^Q^^ 

Ii dt ^ {q + t)2(5g2(2t2 + 1) + l5qt{t^ + 1) + 6t^ + 8t + 1) 6"^* 



ml.^{a,p) = \ > . (106) 



As done in [18] , we present here a qualitative analysis of the ratio of the glueball masses mg 



?n,Q |_. From expressions (104), (105) and (106) we see that, under the hypothesis of Abelian dominance, 

the glueball masses are functions of the threshold parameter a and of the Borel mass //g, encoded in the 
parameter p. According to the SVZ framework, we vary these parameters and we look at the location of 
the glueball masses. The output of our results are shown in Fig.l and Fig. 2. In particular, from Fig.l, one 
can see that, when the parameter a belongs to the interval 1 < a < 1.8, the masses of the three lightest 

states are in qualitative agreement with the available lattice data jl9], i.e. ni^j^j^ < ^i^2++ < ""^o h' ^ 

feature which holds for all values of p, as shown in Fig. 2. As already mentioned, the plots of Fig.l and 
Fig. 2 are precisely of the same kind of those previously obtained in the Landau gauge ^Bj, a fact which 
we interpret as evidence of the gauge independence of the physical spectrum of the theory. Moreover, 
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the existence of the interval 1 < a < 1.8 is regarded as an encouraging consistent confirmation towards a 
more quantitative analysis of the glueballs spectrum in the maximal Abelian gauge [50], along the lines 
presented in [19] . 
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Figure 1: Glueball masses as functions of the threshold parameter a, for p = 5. 
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Figure 2: Glueball masses as functions of p, for a = 1.3 



Let us end this section with the following remarks: 



it turns out that the mass hierarchy for the states 0"*"^ , 2++ , 



the same as that of the Landau gauge, namely r?T,Q++ < m2++ < 



, in the maximal Abelian gauge is 
, and is in agreement with the 

lattice data |49j . We take the present results as evidence of the gauge independence of the physical 
spectrum of the theory, 

it worth to mention that the ordering which we have obtained for the low lying glueball masses has 
also been established in other approaches, such as AdS/QCD and SVZ sum rules. We remind the 
reader to the reference [l2] for a general review on the glueball spectrum. 



the present result can be traced back to two relevant features of the maximal Abelian gauge: i) the 
assumption of the Abelian dominance fTI, '8', "9], ii) the nice fact that the diagonal gluon propagator. 



eq.(50), is of the same kind of that found in the Landau gauge, in both Refined Gribov-Zwanziger 



model [m [221 123] and Landau replica model [20] . 
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6 Conclusion 



In this work we have presented a first quahtative study of the glueball spectrum in the maximal Abelian 
gauge. In particular, the three lightest states, J^^ = 2++, ^ , have been investigated. The output 
of our analysis is that the resulting mass hierarchy is in agreement with both lattice data [l9] and previous 
analytic results obtained in the Landau gauge [18l[19]. As already mentioned, we interpret our result as 
evidence of the gauge independence of the physical spectrum of the theory. 

Concerning the maximal Abelian gauge, we have implemented it within the framework of the replica 
model, already successfully introduced in the Landau gauge [20]. As we have seen, the action of the 
replica model which we have been able to construct in the maximal Abelian gauge enables us to nicely 
reproduce the infrared behavior of the gluon and ghost propagators observed in lattice simulations [17]. 
In addition, the model allows for the introduction of BRST invariant local operators for glueball states 
whose correlation functions display a spectral representation with positive spectral densities. 

We have had also the opportunity to repeatedly state the relevance of the Abelian dominance hypothesis 
OElEj, which is believed to be a key ingredient of the maximal Abelian gauge. We hope that the present 
results might stimulate further investigations on the infrared behavior of the gluon and ghost propagators 
through lattice simulations as well as more quantitative results on the Abelian dominance in the maximal 
Abelian gauge. 
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A Renormalization of the action of the replica model in the maximal 
Abelian gauge 



A.l The starting classical action 



This Appendix is devoted to give a detailed proof of the renormalizability of the action of the replica 
model in the maximal Abelian gauge. Let us begin by reminding the expression of the classical action S 
we start with, namely 

5^ = Slin + Sp + Sfi + S-q + S"^ + S sources ~\~ S extra ■ (10''') 



where Sun , Sp , Sh , S^^ , , S sources ) Sextra are given, respectively, in eqs . (|31 [) , (|34[) , (|38| , (|40[) , (|42[) , (|44[) , (|46 



With the exception of a linear term, irrelevant at the quantum level, the action S displays the following 
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discrete symmetry 

^ y\ 
y ^ y, 

p ^ p, (108) 

where 3^", a = 1,2, y'^'^"-^^ y ^nd p stand for 

yms ^ {A^,Bf,,b,b,c,u;,c,Cj,%,%,L,L 

= {a,a,X,X,p,p,cp,^,h,h,'il),i},r],fi,x,x} , 
p = {m^,M^,fx^,i9'^,v^g,a,Ci,---,C5,ki,...,ke}. (109) 

This symmetry forbids, for example, the occurrence of terms like J df'x a{x)^{x) which are, in principle, 
permitted by the power-counting. 



A. 2 Ward identities 



The action S obeys a large set of Ward identities, which we enlist below: 
• the linearly broken Slavnov- Taylor identity: 



(5c (5L 5lo 6L 5c" (5a)" (5c 6u] SX Sa 

(5S _(5E , (5S -(5S (5S _(5S 5S _(5S 

+p— + + h— + il)-=- + r/— + + + u-= 

d(p op dip oh ox 07] ou ot, 



I 



d'^x (^iV2^'^ X{x) + rr? (f{x) + V'(a;) + x(a;) + u{x)^ . (110) 



the equations of motion of the diagonal Lagrange multipliers and of the antighost equations: 

the equations of motion of the auxiliary fields: 

§ = -{a-iV2^'), ^ = -A; (113) 

^ = -(p-m^), - = -^; (114) 
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" Su 

• the linearly broken parametric equations: 



({-.^), g = -«. (117) 



as 

gs 
as 
as 

the linearly broken integrated Ward identity: 



xa{x); (118) 

(119) 

Jd'^xh{x)- (120) 
y" ct^xf]{x) ; (121) 
Jd^x^x). (122) 



TS = d^x{—-a—-ic''—-iuj''-=- + 2—— + 2—-y- 
J \d(p Sip Sb"- Sb'^ Su 5L du 6L 

= J d^x {<f - ai>) . (123) 

• the diagonal local ghost equations: 

= -a„(a„c+!!^)+9£'*(!i;4-'^°<:'); (i^i) 

= -d^{d^u + ^^) + ge''\niB'^-L''u'>). (125) 

• the 5^(2, SR) Ward identities: 

P(S = / d^ic c»— -z— — + 2— =0; 126 

[a f (5S (5S (5S(5S\ , 

PS) = d^xloj"- i-^-^ + 2-^-=] =0. (127) 
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the linearly broken U{1) residual identities, typical of the maximal Abelian gauge: 

(5S 



a 



-id% + fi^d^A^; (128) 



ME) . a„j|,8„(B„|)...-|:./| 



= -id^ + fx^d^Bf,; (129) 

e ^ {5^,6^cD^^^^7;l,Z'^}. (iso) 

A. 3 Characterization of the most general counterterm 

In order to characterize the most general local invariant counterterm compatible with the whole set of 
Ward identities we follow the procedure of the Algebraic Renormalization [^Slj. We perturb thus the 
action S by adding an arbitrary local integrated polynomial E'^"""* in the fields and external sources, of 
dimension bounded by four, and we require that the perturbed action, S + eS'^"""*, obeys, to the first 
order in the expansion parameter e, the same set of Ward identities fulfilled by S. This gives rise to the 
following constraints for J]™""*: 

5eS^°""* = 0; (131) 

TeS'=°""* = 0; (133) 

Gc S™™* = , Gu; = ; (134) 

= , S'^"""* = ; (135) 

WS'=°""* = 0, WS'=°™* = 0; (136) 
where stands for the linearized nilpotent Slavnov- Taylor operator 

f . f 6^ 5 6^ 6 6^ S 61: 6 6^ 6 6^ 6 6^ 6 61: 6 



with 



+ ^f^JTT = ' S^"""* = , (132) 



J K6A-6n- 6^1 6 A- 6B-6n- 6n-6B- 6A^6nf, 6%6A^ 6B^6Vl^ 6%6B^ 
<5S 5 6J: 6 61: 6 61: 6 61: 6 61: 6 61: 6 61: 6 6 6 

^J^Ji^^ Ji^I^^ I^Jl^^ Jl^J^^ J^Jl^ JlTc^ J^Jl^ JlJ^^'^ 6^^^ 6o^ 

,6 - 6 6 - 6 6 6 6 - 6 6 6 6 6 \ , 
+ib- + lb— + 'yj^ + >^j^ + + ^j-- + h— + + + Xjz + + u-=] , 137 
oc duj OA OCT dip dp dip oh dx 6r} du 6t,J 

S^S^ = , (138) 

and Ts , , V^. are given by 

f ,A ( 6 6 „ 5 „ 5 61: 6 61: 6 61: 6 61: 6 \ , , 
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(Tx c" 



i—— + 2—-^ + 2 



d^x - i^^- + + 2^^ 



(140) 
(141) 



From the equations of motion of the Lagrange multipHers, eqs.( 111 ),( 112 ), and from those of the auxihary 



fields, eqs.(113)-(117), it fohows immediately that the counterterm S'^"" must be independent from the 
set of fields 

{h,h,a,\p,^,fi,x,h,i},i,u] . (142) 

Similarly, from the linearly broken parametric equations, eqs.( 122| ), it follows that S™""* is independent 
from the mass parameters i?^, m^, M^, /x^ and v"^. Due to the nilpotency of the linearized Slavnov- 



Taylor operator eq.(138), the characterization of the most general allowed counterterm corresponds 



to identifying the cohomology [31j of the operator Sy, in the class of the integrated local polynomials of 



dimension four and zero ghost number, subject to the remaining constraints, eqs.(132)-(136). After a 



lengthy algebraic analysis, it turns out that the most general invariant counterterm can be written as 



^count 



ao Sym[A] + ao Sym[B] + cS^A 



(-1) 



(143) 



with '^^ being given by 



MM 



+\ (02 + 03 + a as) V^(A;1^^ + B^Bl) - 0(03 + a4 - ag) (/^(c'^c'^ + cD'^^^'^) 
+^ (as + ay A + ag x)(^;iA^". + 5^5^) + (ag V + ^8 A + aio x)(c"c'^ + c^'^w") 

I ab f —a b I —a b\ \ Cl \ , Cs / 7 . C4 . Cs /- 

+aii s-e ^(c c + w ) + — ai2 Act + y ai3 + y vh+ — ais X?? + y < 



+/ci ai7 Ap + k2 ai8 A?? + A;3 aig A/i + /C4 020 '^'H + ^5 a2i + /cg 022 



(144) 



and Qn, with n = 0, . . . , 22, being arbitrary dimensionless coefficients. 



It remains now to show that the counterterm (143) can be reabsorbed into the starting classical ac- 



tion E through a suitable redefinition of the fields, sources and parameters. This will be the task of the 
next subsection. 



A. 4 Determining the renormalization factors 



In order to show that the counterterm J]'^" , eq.(143), can be reabsorbed through a redefinition of the 



fields, sources and parameters, we shall employ the following compact notation. All fields will be denoted 
by {/}, all sources by {J} and all parameters by {p}. As usual, the renormalization factors are introduced 
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by 



fOff 

JO 

/'diag 




j; 



off 



j; 



diag 



= zfrff, 

yV^ fdiag 

- ^f J 

= zjrff, 

= Zj J'^'"-3 , 



Pq = Zpp, 



(145) 



where the quantities labeled by "0" stand for the bare, or unrenormalized, quantities, while those with no 
label are the renormalized ones. Stating that the counterterm S^o""* can be reabsorbed into the starting 
action E, amounts to prove that 



S[/o, Jo,Po] = J,p]+s J,p] + 0{e^) 
Due to the mirror symmetry, it follows that 



(146) 



Za = Zu , Za = Zu 

Zp = Z,T, , Zr = Z/7, , 



Zi) , Zh 



Zn 



Zn 



-'n ) 



Zc- 
Zl 



Z,, 



z.-- 

Zl 



Zu 



(147) 



For what it concerns the renormalization factors of the linearizing auxiliary fields and of the mass param- 
eters, wc have to properly take into account that these variables can mix each other, due to the fact that 
they possess the same quantum numbers, i.e. the same dimensions, the same color structure , etc. As a 
consequence, the corresponding renormalization factors turns out to be expressed in terms of matrices, 
namely 



(148) 



(po\ 






























ho 




h 




ho 




h 




V'o 


= Z3 






^0 


= Z4 






= Zi 


a 






= Z2 


a 




Ao 


A 




Ao 


A 






\v) 




Km) 








\Xo) 








\XoJ 







(ml\ 






Ml 


= Z5 











(149) 



All the renormalization factors, even the Z matrices, coincide with unity at zeroth order in the pertur- 
bative expansion. 



Z 

z 



l + 0(e), 
1 + 0{e) . 



(150) 
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After direct inspection, it turns out that condition (146) can be consistently fulfilled with the Z's factors 
determined in the following way: 



and 



z 



1/2 



Za 

Zki 

Zk4, 



1 -e 



ao 



1 + e 



02 _ 03 
2 2 

l + .(^ + ^ 
'22 



ao 



a2 + Oil 



l + e 

l + e 
l + e 
l + e 
l + e 
l + e 
l + e 



l + e [an — 02 — as — 005 
l + e (ai8 - 2ao) , 
l + e {aig - ao + 02 - ae - 

l + e ^a2o — a2 — 03 — aas - 



(ao — 2a2 + a^) , 
(ai2 + 2ao) , 

(ai3 + 2ao — 2a3 — 2aa5) 
(ai4 + 2a2 — 2a6) , 
(ai5 - 2ao) , 

(ai6 - ao + 2a2 - 2aii) , 



07 



l + e 



021 + ao - (1 



l + e ^a22 — ao + 02 



ae - 



2C2 

h 
2C3 , 

as + — 04 



ag 



2C2 , 

a - — as 
K5 



1 + 1— 06 



2C_3 
fee 



aio 



l + e 



/— ao + a2 + as + aa^ 
—a{as + 04 — ae) 



V 









as 


a? 


ag 


a2 + ae 


as 


aio 





ao 











ao 



(151) 



(152) 



l + e 



ao — a2 — as — aas 
-as 
-aj 
-ag 



a(as + a4 — ae) 
02 - ae 
-as 
-aio 















-ao 








-ao 



(153) 
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1 + e 



2 v^-u - 02 — as) + "05 
—0(03 + 04 — ae) 





05 

2(00 - 3a2 - as) + ae 





07 
as 

2a7 — 00 — 02 — 03 




aio 


4(ao + 02 + 03) 



(154) 



/^(ao 



1 + e 



a2 - 03) - 005 
-05 
-07 
-09 



a(a3 + 04 — oe) 
-5(00 - 802 - 03) - Oe 



-oio 






-2o7 + 00 + 02 + 03 











5(00 + 02 + 03)^ 

(155) 



1 + e 



00 + 02 + 03 + aos 


05 


i\/2 07 


09 


— a(o3 + 04 — Oe) 


—02 + 06 


i-v/2 og 


OlO 








Oo 














Oq 



V 

The remaining renormalizations factors are not independent, being given by 



Z 



1/2 



^1/2 ^ 02 ^ 



^1/2 ^ ^-1/2 



■A/2 



y 71/2^1/2 



1/2 



7 _ 7-1 7-1/2 7-1/2 



z^ = z]!\ 



= z-^^z-'i'~z-'i' 



^1/2 ^ ^-1/2 



1/2 



7I/2 



7-171/2 7I/2 
9 -^c -^^ 



Zi^ = z-^^z-^ , 

1/2 _ ^-1/2 ^-1/2 
u — ^g^c 



(156) 



This proves that the action S of the rephca model in the maximal Abelian gauge is renormalizable. 



(157) 
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